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It is of great interest whether the recent KamLAND bound on the flux of electron antineutrinos 
from the Sun constrains the electromagnetic properties of the neutrino. We examine the efficiency 
of the electron antineutrino production in the solar magnetic fields, assuming the neutrinos are 
Majorana particles with a relatively large transition moment. We consider fields both in the radiative 
and convective zones of the Sun, with physically plausible strengths, and take into account the 
recently established values of the oscillation parameters. Our analysis shows that the production rate 
in question is presently unobservable. In the radiative zone, it is suppressed by the large measured 
value of the flavor mixing angle which eliminates the resonant level crossing. A corresponding general 
resonance condition, valid for large as well as small values of the mixing angle, is derived. Likewise, 
in the convective zone, the strength of the small-scale magnetic field is likely to be insufficient. Thus, 
no useful bound on the neutrino transition moment can be derived from the published KamLAND 
bound. KamLAND may be, however, on the edge of probing an "optimistic" scenario, making 
further improvements of its sensitivity desirable. 

PACS numbers: 13.15.+g, 26.65.+t, 14.60.Pq 



I. INTRODUCTION 

One of the most important current tasks in neutrino 
physics is to understand the full implications of the so- 
lar HSU EMM 11 IS 113 , atmospheric HI G2, 
and reactor |13l Il4| neutrino data that became available 
in the last five years. While much of the recent effort 
has been focused on establishing the values of the neu- 
trino masses and mixings, there has also been a growing 
realization that the same data may contain valuable in- 
formation about other neutrino properties, such as its 
interactions with matter [TH ITfi ITtI Eft pfl l2ll pp 
or with electromagnetic fields |23l 1241 l2a 12a. 1271 l28l |29| . 
Any measured deviation of these properties from their 
Standard Model (SM) values would have profound theo- 
retical implications. 

As part of this program, we wish to reexamine here the 
sensitivity of the recent solar neutrino measurements to 
the neutrino electromagnetic properties. 

The idea that the neutrino may possess a relatively 
large magnetic (transition) moment has been the sub- 
ject of intense theoretical speculations for many years. It 
was originally put forth in the early 1970's, as a possible 
explanation to the solar neutrino deficit 31]: the solar 
neutrinos, it was argued, were "disappearing" as a re- 
sult of the neutrino spin precession in the solar magnetic 
field. Remarkably, this ex pla nation - much improved 
with time |32l 133 . 134 1351 1361 137| - remained viable for the 
next three decades|69^Qnlv in 2002 did the KamLAND 
measurement of the reactor antineutrino disappearance 
|13| conclusively rule it out as the dominant mechanism 
of the solar neutrino conversion. The possibility of the 
spin precession happening at a subdominant level, how- 



*Electronic address: friedland@lanl.gov 



ever, remains of great interest to this day, as a probe of 
the neutrino electromagnetic properties and, at the same 
time, of the magnetic fields in the solar interior. 

For definiteness, here we will focus on the case of Majo- 
rana neutrinos. The experimental sensitivity to magnetic 
spin precession in this scenario is especially high, as has 
been long recognized |38| : upon the spin flip, a Majorana 
neutrino turns into an antineutrino and - thanks to flavor 
oscillations - can then be detected with some probability 
as an electron antineutrino, v e . The v e flux from the Sun, 
in turn, is strongly constrained by Super-Kamiokande |(| , 
SNO |a, and especially KamLAND 8] which can detect 
v e — ► v e conversion of the solar 8 B neutrinos at the level 
of a few hundredths of a percent. 

To understand if the KamLAND bound places a non- 
trivial constraint on the neutrino magnetic moment, we 
need to compare it to the maximum possible v e flux from 
the Sun. This must be done taking into account (i) the 
recently measured values of the oscillation parameters, 

(ii) the bounds on the solar magnetic field strength and 

(iii) the existing direct bounds on the neutrino magnetic 
moment. We also need to consider that magnetic spin 
flip could occur either in the convective zone (CZ) of the 
Sun, which comprises the outer 30% by radius, or deeper 
in the radiative zone (RZ). In the first case, the neutrino 
traverses a region of turbulent magnetic field; in the sec- 
ond case, the field, if present, is smooth and frozen into 
stationary plasma. The physics of the neutrino evolution 
is quite different in the two cases and requires separate 
treatments. 

Both the cz mmm and the rz mm us cases 

have been investigated in the literature and in both cases 
it has been concluded that it was a priori possible to have 
a signal greater than the current KamLAND sensitivity 
reach. In the case of the CZ, even a bound on the neu- 
trino magnetic moment was claimed. 

In our analysis, we reach different conclusions. In the 
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case of the RZ, we find that the usual treatment, based 
on the v e — ► resonance, works only for small flavor 
mixing. The large measured value of the solar neutrino 
mixing angle leads to a qualitatively new effect: the res- 
onant level crossing disappears and no measurable mag- 
netic spin flip can take place. A new analytical result, 
extending the classical "magnetic resonance" condition 
|35ll36j | to large flavor mixing, is derived here for the first 
time. 

In the case of the CZ, we consider two qualitatively 
different models of the magnetic field. One of these mod- 
els, that of the uniform Kolmogorov turbulence, follows 
[28l l29j| and our basic treatment of the neutrino evolu- 
tion in this case agrees with that work. We, however, 
take what we regard as more physically justifiable mag- 
netic field parameters. For the second model, we take 
a scenario of isolated flux ropes as suggested by the re- 
cent solar physics research. This scenario, to the best 
of our knowledge, is considered in the neutrino literature 
for the first time. The net result is that in both cases we 
estimate the maximum v e flux to be below the present 
KamLAND bound. 

We begin by giving the necessary background informa- 
tion, including the form of the interaction (Sect. Ill A[) . 
the known bounds on the neutrino magnetic moment 
fSect. ITT*B)l . the general properties of the solar magnetic 
fields (Sect. and the form of the neutrino oscillation 
Hamiltonian fSect. ITT"D|) . The case of the RZ magnetic 
field is treated in Sect. IIIII after reviewing the conven- 
tional approach fSect. lTTl A|) . we argue that resonant an- 
tineutrino production is suppressed for large flavor mix- 
ing (Sect. ITIlBI . derive a corresponding general analyt- 
ical condition (Sect. ITTl Cf) . and finally perform a direct 
numerical scan of the parameter space to confirm our an- 
alytical conclusions (Sect. ITiTD|| . This is followed by the 
analysis of the CZ fields in Sect. II VI where the two mod- 
els of the magnetic field mentioned above are considered 
(in Sects. IIV Al and IIV B|) . The results are summarized 
in Sect. [V] Some additional technical details about the 
CZ neutrino evolution can be found in the Appendix. 



II. BACKGROUND 

A. Interaction 

The electromagnetic moment interaction of the Majo- 
rana neutrinos is given by the following dimension five 
operator 

CeM = -^Hab^a^^bF^, (1) 

where and 'J'b are Majorana spinors of flavors a and 
b in the 4-component notation, and E'*" = £[7'*, 7 I/ ]/2, 
with /i and v being the Lorentz indices. 

It proves useful to rewrite this Lagrangian in the 2- 
component Weyl spinor notation. This notation, besides 



being simply more natural for ultrarelativistic fermions, 
brings considerable physical insight. 

Writing the 4-component Majorana spinor 4" in the 
Weyl basis as *Sf = (Xa,X a ) T an( i introducing 2x2 ma- 
trices (o-n a = o^**" ")/2, with = 

(1, a), er Q/3 v = (1 , — (?) (here a are the usual Pauli matri- 
ces), we rewrite Eq. as 

£em = -^ ab (x a )a(an a f3 (xp)bF^ + h.c. (2) 

This form makes the following important physical prop- 
erties of the interaction manifest: 

• The interaction couples the neutrino and antineu- 
trino states. Notice in Eq. the spinors are either 
both "undotted" or both "dotted". Each "undot- 
ted" spinor destroys a neutrino state or creates an 
antineutrino state. 

• No flavor-diagonal terms. Since spinors anticom- 
mute, the a = b terms in Eq. vanish identically, 
and hence Majorana neutrinos cannot have flavor- 
diagonal magnetic moments [sil Isi^ . This can be 
easily seen, for example, using the definition of the 
(0-^)0? matrix given above and the identity 

x<&r$ = Tp<r u <^x (3) 

The transition (a ^ b) moments, however - which 
lead to a simultaneous spin and flavor change, i.e., 
v e —> T - are perfectly allowed. 

• The interaction with the z (longitudinal) compo- 
nent of the magnetic field vanishes identically, 
while the interaction with the transverse compo- 
nents B± = (B x ,B y ) is nonzero. One of the two 
spinors in Eq. (j2J has a raised index. Denoting it 
by x = X a — eQ ^X/3? we write Eq. J5J in the three- 
dimensional notation as 

Cem = -Vab{x)ad{v)b{B - iE) + h.c. (4) 

Plugging in Xo « Xi « (0 I) T - which neglects 
terms of order m/E u - we find that the interaction 
picks the 12 component of the matrix a, i.e., 1 for 
a x , —i for o~ y , and for a z . 

The same conclusion can be reached by boosting 
the magnetic fields to the neutrino rest frame. 

• The real part of fi a b corresponds to the magnetic 
dipole moment and the imaginary part to the elec- 
tric dipole moment. An ultrarelativistic neutrino 
precesses in an external magnetic field if either the 
magnetic or electric dipole moments are non-zero 
[33l l34| . Following the convention, we will hence- 
forth use the term "magnetic moment" to refer to 
both possibilities. 

All these properties are well-known and are derived 
here for completeness. 
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B. Size of the transition moment 

The laboratory bounds on the neutrino magnetic 
(transition) moment come from measuring the cross sec- 
tions of ve~ or ve~ scattering in nearly forward direction. 
The recent bound for the interaction involving the elec- 
tron antineutrino is /i e < 1.0 x 10~ 10 /j,b at the 90% con- 
fidence level j^jji where [ib = e/(2m e ) is the Bohr mag- 
neton (m e is the electron mass, e is its charge). Stronger 
bounds, // < 3 x 10~ 12 /is, exist from astrophysical con- 
siderations, particularly from the study of red giant pop- 
ulations in globular clusters ^l|. Larger values of the 
transition moment would provide an additional cooling 
mechanism and change the red giant core mass at he- 
lium flash beyond what is observationally allowed. 

Even allowing for new inputs (such as updated dis- 
tances to globular clusters) or refinements in the stellar 
models, it appears extremely unlikely that the red giant 
bound could be weakened by more than a factor of a few. 
Correspondingly, as a very conservative bound, we will 
take /i < 10~ 11 fiB in the subsequent analysis. 

From the theoretical point of view, the neutrino does 
possess a nonzero transition moment, already in the SM. 
It is induced by loop effects: as the neutrino briefly dis- 
sociates into virtual charged particles (W and a charged 
lepton I), those constituent particles couple to the pho- 
ton. Because the operator in Eq. J2J) requires a helicity 
flip, and because in the SM the W boson couples only 
to left-handed fields, the mass insertion has to be put on 
the neutrino line. The resulting moment is proportional 
to the neutrino mass and hence is highly suppressed, 
/i„ ~ eGF"m v ~ 10 _19 /is(TO t/ /ey). In possible exten- 
sions of the SM, however, this need not be the case, and 
the effect can be proportional to the mass of the charged 
lepton running in the loop. For example, in the left- right 
symmetric model one obtains \i v ~ eGpnii sin 2rj |42fl . 
where r\ is the left-right mixing parameter in the model. 



C. Solar magnetic fields 

1. Convective Zone 

The CZ of the Sun is known to host significant mag- 
netic fields, which manifest themselves in a variety of ef- 
fects on the solar surface (sunspots, flares, prominences, 
etc). These fields are created and destroyed during each 
solar cycle by a combination of convective motions and 
differential rotation. Detailed understanding of the solar 
cycle is still an active subject of research. Nevertheless, 
the following general features of the CZ fields are well 
established: 

• The magnetic field observed in sunspots is several 
kG, exceeding the average surface field strength by 
three orders of magnitude. Sunspots usually come 
in pairs of opposite polarity and are thought to 
be manifestations of large-scale magnetic structures 



residing inside the CZ. 

• The total toroidal flux in the CZ at sunspot max- 
imum can be estimated from the total flux that 
emerges on the surface during the solar cycle, which 
is around 2 x 10 25 Mx ^3|- Since the same flux tube 
may be emerging more than once, 10 24 Mx has been 
argued to be the likely value [43j . Averaged over 
the CZ, this gives a field of several kG. 

• The turbulent equipartition value for the magnetic 
field at the base of the CZ can be expressed in 
terms of the density p, the solar luminosity Lq, 
and the distance to the center of the Sun r: B eq ~ 

• Modern models of the CZ fields argue that fields on 
the order of 100 kG may exist in the CZ. For the 
present discussion, it is important to realize that 
such strong fields, if they exist, can occupy only a 
small fraction of the CZ volume, such as the thin 
shear layer near the base of the CZ ('ta chocline'), 
or be localized in isolated flux tubes JSj . To have 
even half of the CZ filled with such strong fields 
would be contrary to the flux arguments mentioned 
above, as well as basic energetics considerations. 

For an in-depth review of the CZ fields and further 
references, the reader is urged to consult an excellent 
review |44j . For the status of the present-day simulation 
efforts see, e.g., |45|. 

2. Radiative Zone 

Unlike the CZ, the RZ has no internal motions to pro- 
vide an active mechanism to generate magnetic fields. 
On the other hand, the conductivity of the RZ plasma 
is very high, which opens up the possibility that the RZ 
could support large-scale magnetic fields left over from 
the early stages of the solar evolution. Let us consider 
both of these points in turn. 

The evidence for the lack of any large-scale mixing mo- 
tions in the RZ comes from several sources. First, helio- 
seismological studies show that, unlike the CZ, the RZ 
rotates as a solid body. (This fact by itself, inciden- 
tally, provides a hint that at least a weak poloidal field is 
present in the RZ [46j). Second, solar model calculations 
show that the stratification (entropy gradient) in the RZ 
interior is quite strong, precluding large radial motions. 
Finally, a very impressive direct confirmation of the lack 
of radial mixing is provided by the observed Beryllium 
abundance at the solar surface [47|. Be is destroyed at 
T = 3 x 10 6 K, which is the temperature at r = O.61i?0, 
Rq being the solar radius, 6.96 x 10 5 km. Since Beryl- 
lium observed at the solar surface is not noticeably de- 
pleted relative to its primordial abundance, the material 
at r ~ O.7i?0 (on the bottom of the CZ) and the material 
at r ~ O.6i?0 have not been mixed during the lifetime of 
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the Sun (except for a brief period during the pre-main 
sequence convection stage). 

While a large-scale magnetic field in the RZ is not be- 
ing presently generated, it is also not destroyed by Ohmic 
decay over the solar lifetime, 4.7 billion years. The Ohmic 
decay time for large-scale field configurations is in the bil- 
lions of years 0] , in fact for the lowest toroidal mode it 
is about 24 billion years |49| . Over time, a primordial RZ 
field has been invoked by many authors to explain a va- 
riety of effects: from the solar neutrino problem [5(3, [ED 
to the so-called Princeton solar oblateness measurements 

What is important for the present analysis is that the 
field in the RZ, if it exists, must be (i) smooth and (ii) 
bounded in strength ^j. Smoothness follows from the 
fact that, because of Ohmic decay, any small-scale mag- 
netic field (I < Rq/10 — Rq/20) features would have de- 
cayed away over the lifetime of the Sun. This ensures the 
adiabatic character of the neutrino evolution throughout 
the bulk of the CZ. The bounds on the field strength 
come from several independent arguments, such as the 
measurements of the solar oblateness and the stability 
analysis of field configurations. It follows that the ampli- 
tude of the large scale toroidal field cannot exceed 5 — 7 
MG |j[J . Combined with the bounds on the size of 
the neutrino transition moment, this limits the value of 
the magnetic term (see below) in the neutrino oscillation 
Hamiltonian. 



D. Oscillation Hamiltonian 

Since flavor oscillations change the neutrino (and an- 
tineutrino) flavor and the magnetic interactions convert 
neutrinos into antineutrinos, the evolution of the neu- 
trino state is, in general, governed by a 6 x 6 Hamiltonian 
matrix, with the neutrinos and antineutrinos of all three 
flavors. For the purpose of studying the spin flip in the 
RZ, the problem can be reduced to a 4 x 4 case [jlj. I n 
the basis (y e , v^v e ^v^)^ the Hamiltonian is [sfj 



H = 



H v (B x 

¥ %By)M 



- iBy)M^ 



(•5) 



where B x y are the transverse components of the mag- 
netic field and the 2x2 submatrices are given by 



M 

Hr, 



-^e M 

-A cos 29 + A e 
A sin 29 

-A cos 29 - A e 
A sin 29 



A sin 29 
A cos 29 + Af, 

A sin 29 
A cos 29 - A b 



(0) 
(7) 
(8) 



Here A = Am 2 /4E l/ , A e = y/2G F ( 

n e — n n /2) and 

A^ = \/2G p(—n n /2). Am 2 is the neutrino mass-squared 
splitting, E v is its energy, and n e and n n are the electron 
and neutron number densities. 



The physical reason for the 6^4 reduction is that the 
large atmospheric neutrino splitting effectively decouples 
two of the six states in the case of the smooth magnetic 
field. It should be noted, however, that the "noisy" field 
in the CZ enables transitions between states separated 
by the large splitting. Consequently, in the CZ one may 
need to consider the full six- neutrino problem |54j |. 



III. RADIATIVE ZONE 

A. Conventional approach 

Historically, in discussions of spin-flavor precession 
(SFP) in the Sun, special attention has been given to the 
"SFP resonance" region [3^, |3(|, defined as the region 
where two of the diagonal terms in the Hamiltonian 10- 
|SJ, namely those corresponding to the ee and JipL entries, 
become equal. The corresponding resonance condition is 



V2G F (n e -n n ) = — ^- cos 29. 



2E V 



(9) 



Compared to the traditional "flavor resonance" condi- 
tion, defined as the equality of the ee and fifj, diagonal 
entries in H5I8|) ~). 



V2G 



F n e 



2E„ 



■cos 20, 



(10) 



condition JjJJl is satisfied at a somewhat higher density in 
the Sun. 

When the resonant SFP mechanism was introduced 
back in the late 1980's, the values of the oscillation pa- 
rameters were essentially unknown. In contrast, nowa- 
days the solar and KamLAND experiments have told us 
that Am 2 - 8 x 10~ 5 eV 2 and tan 2 9 ~ 0.4 pf. Two im- 
plications of this have been recognized: (i) if the process 
of SFP operates in the Sun, it must be combined with the 
"ordinary" flavor oscillations [H HI |H H3 , and ^ the 
SFP resonance condition, as given in Eq. 0, is satisfied 
for 8 B neutrinos in the RZ. 

The second point is illustrated in Fig. ^ where the 
four diagonal elements of the Hamiltonian arc plotted as 
a function of the distance from the solar center. Here and 
later in this paper, the BP04 solar model J5f| is used. 

In most of the recent studies of the RZ spin-flip [24II27I 
l30j . the two processes, flavor oscillations and spin-flavor 
precession, are treated as occurring in spatially separated 
regions. The resulting probability of v e — + v e conversion 
is then given as a product 



P{y e -> D e ) = P(u e 



.) x P(i 



(11) 



A notable exception is Ref . [23] , where it is noted that 
the evolution in the case of the RZ does not decouple into 
two separate problems. The numerical conclusions pre- 
sented in Fig. 2 of that paper, however, are quite different 
from ours, as discussed later. 
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FIG. 1: The diagonal elements of the Hamiltonian 15181 for 
Am 2 = 8 x 10" 5 eV 2 , E u = 9 MeV. The solid/dashed lines 
correspond to the neutrino/antineutrino entries. 



B. Large mixing and the disappearance of the 
resonance 

The above treatment relies on the notion of resonance, 
as given by Eq. @. Fundamentally, this notion comes 
from the analysis of a two-level system, and even there 
it applies only in the limit of small mixing [56L l57j . We 
wish to examine whether it can be applied to the four- 
level system at hand. 

Let us first briefly review the ph ysic s behind the stan- 
dard small-angle MSW effect [53, |5!j ■ Suppose we are 
given a two-level neutrino system with a Hamiltonian 



-D d 
d D 



(12) 



If D ^> d, flavor oscillations are suppressed, as observed 
already in |58j . Indeed, the mixing angle is given by 
tan 20 — d/D and the oscillation amplitude is propor- 
tional to sin 2 20. The physical idea behind the resonance 
is that this suppression is lifted if the diagonal splitting 
becomes small in some region along the neutrino trajec- 
tory, as a result of the crossing of the levels. Suppose 
|-D(/)| < \d(l)\ over a distance 61. Then, inside this in- 
terval, d becomes the dominant part of the Hamiltonian 
and hence drives oscillations. If 51 is sufficiently large, it 
may be possible to obtain large flavor conversion |72j . 



With this in mind, let us return to the Hamiltonian 
in Eqs. 1|5I8|) . Now, the electron neutrino state, in addi- 
tion to being coupled to the state via the magnetic 
moment interaction, is also coupled to state via the 
flavor oscillation term. A simple, but crucial observation 
is the following: for large values of the mixing angle, even 
if the diagonal ee and ftp, terms are equal, the Hamilto- 
nian is still not dominated by the off-diagonal magnetic 
terms /iB, but by the off-diagonal flavor oscillation terms 
A sin 29. Thus, there is no reason to expect that the 
condition in Eq. © should automatically lead to large 
v„ — > v u conversion. 



This complication is, however, easily removed by going 
to a basis in which the large off-diagonal flavor oscillation 
terms are absent. Consequently, we diagonalize H v and 
-ffp given in Eqs. J7J and (|SJ by performing separate flavor 
rotations on neutrinos and antincutrinos. Wc find 



H 



I Ami 



Me/x-B sin(0 m 

\flenB COs(0 m 





A m 2 

Ojn) -He^,BcOs(9 r , 

6m) V-e^B sin(# m 



li eiJi Bsm(9 m 
-fx efi Bcos(8 r 

A m i 





- 9m) He^B COs(6» m 

-0 m ) n ef *Bcos{8 m 


A m2 



(13) 



r 



In Eq. (pil, A m i, A m2 , A m i and A m2 are the eigen- 
values of the oscillation Hamiltonian (15181 for vanishing 



magnetic field, 



± 



A cos 29 



t _ A e + A n 

L^mi — ^ ^ 



A,? A* 



A cos 29 



r 



A. P A. r , 



A 2 sin 2 29 



A 2 sin 2 29 



(14) 
(15) 



while 9 m and 9 m are the rotation angles, The barred quantities refer to antineutrinos. 

tan26» = A sin 20/ [A cos 20 - (A e - A n )/2], (16) 
tan20 = A sin 20/ [A cos 20+ (A e - A n )/2]. (17) 
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Let us first discuss the general properties of these equa- 
tions. When the matter potentials are negligible, A ml 
and A m i both approach —A, while A m 2 and A m 2 both 
approach A. Likewise, in the same limit the rotation 
angles 9 m and 9 m both approach the vacuum value 9. 
Notice that in vacuum the structure of the off-diagonal 
magnetic interactions in the mass basis is the same as 
in the flavor basis. This structure is that of an antisym- 
metric tensor, Eq. ©, which is invariant under an SU(2) 
rotation applied to both indices. More generally, this is 
a manifestation of the fact that the operator in J2J) can- 
not convert a Majorana neutrino state in vacuum into its 
own antiparticlc, whether this state is a mass eigenstate 
or a flavor eigenstate. In sufficiently dense matter, how- 
ever, 9 m and 9 m are, in general, different, because the 
matter terms for neutrinos and antincutrinos have oppo- 
site signs. Thus, in matter, direct transitions between all 
four states become allowed 73] . 

For the energies at which KamLAND is sensitive to 
solar antineutrinos, E v > 8.3 MeV, the magnitude of 2A 
is of order - 8 x l(T 5 /(2 x 10 7 ) - 4 x 1CT 12 eV. By 
comparison, the value of the off-diagonal magnetic terms 
is at most of order lO^V-B x 5 MG ~3x 10~ 13 eV, or 
at least an order of magnitude smaller. This makes the 
analysis in the mass basis quite straightforward. Clearly, 
transitions between any two states with diagonal splitting 
of order 10~ 12 — 10 -11 eV will be suppressed. Conversions 
between a pair of states with smaller diagonal splitting 
may take place, provided that 9 m and 9 m are sufficiently 
different, i.e., the matter potential is sufficiently large. 

In the top panel of Fig. |2 we show the diagonal ele- 
ments of the Hamiltonian Eq. (|13|l for the best-fit values 
of the oscillation parameters. Obviously, their behavior 
is quite different from the diagonal elements in the flavor 
basis, plotted in Fig. ^ I n particular, in Fig. [5] there is 
no level crossing at R ~ O.18i?0. The mass eigenvalues 
do not intersect anywhere inside the Sun. 

Notice that the top two levels in Fig. [21 get close to 
each other for r > 0.3 — 0.4i? Q . The difference 9 m — 9 m 
there is small, but non-zero, and a priori could drive a 
transition between these two levels. The possibility of 
this transition requires a more detailed study. 

The answer comes from considering the adiabaticity of 
the neutrino evolution. Consider an idealized problem: 
a neutrino mass eigenstate moving through a magnetic 
field which is adiabatically turned on and then turned off. 
The neutrino will be partially converted into the corre- 
sponding antincutrino state while inside the field region, 
but upon exiting the region will revert back to being fully 
a neutrino. 

This idealized description applies in our case: (i) At the 
production point in the core, for the best-fit KamLAND 
mass splitting Am 2 ~ 8x 10 -5 eV 2 , the diagonal splitting 
dominates, suppressing the effects of the magnetic field. 
The neutrino is produced predominately in the heavy 
eigenstate (the top line in the upper panel of Fig. EJ) . (ii) 
In the region r ^ 0.2 — 0.5R@, the magnetic term drives 
partial v <-> v conversion (transitions between the top 
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FIG. 2: The diagonal elements of the Hamiltonian (11311 (top 
panel) and the mixing angles 8 m and 9 m (bottom panel) for 
Am 2 = 8 x 1(T 5 eV 2 , tan 2 9 = 0.4, E v = 9 MeV. The 
solid/dashed lines correspond to the neutrino/antineutrino 
levels. 



two states in the top panel of Fig. 01. (iii) As the effects 
of the magnetic term are decreased (mainly due to the 
alignment of 9 m and 9 m at smaller densities), the v v 
conversion is largely undone. The adiabaticity through 
the bulk of the RZ zone is ensured by the physics of 
the RZ fields: as mentioned in Sec. Ill C 21 any structures 
smaller than i? Q /10 — i? Q /20 must have decayed away 
over the lifetime of the Sun. 

This behavior is illustrated in Fig. [21 which shows that 
the electron antincutrino fraction P indeed first grows, 
but then decreases to unobservable levels, as the effect of 
the magnetic term decreases in the outer part of the RZ. 
The residual nonzero value of P is due to the violation 
of adiabaticity near the top of the RZ. 

In summary, except in the limit of small mixing angle, 
Eq. @ is not the right condition for describing resonant 
neutrino-antineutrino conversion. A correct physical cri- 
terion should be based on comparing the diagonal ele- 
ments in the mass basis, not in the flavor basis. This 
"mass basis level-crossing" would be indeed realized in 
the Sun, had the value of the mixing angle been small [74|. 
For the large measured value of the mixing angle, how- 
ever, the level-crossing is simply absent and the efficiency 
of the neutrino-antineutrino conversion is very low. 
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Q_ 




FIG. 3: A typical case of the neutrino evolution, for the value 
of the RZ field near the limit. Shown is the probability of 
finding a neutrino in the z? e state as a function of the distance 
from the solar center (in units of the solar radius), for Am 2 = 
8 x 1(T 5 eV 2 , E v = 10 MeV, 6 = vr/6, B max = 5 MG. 



C. Large mixing and level crossing: analytical 
treatment 

Let us now derive the general analytical condition for 
having resonant neutrino-antineutrino conversion. Refer- 
ring to Fig. |3 we see that, a priori, either the two light 
levels or the two heavy levels could intersect. The corre- 
sponding analytical conditions for these two possibilities 
are 

A m 2 = A m 2, (18) 

and 

A m i = A ml , (19) 

where the indices "1" and "2" refer to the light and 
heavy states respectively (the "-" and "+" signs in 

Eqs. const ). 

Written out explicitly, the above two conditions read 



aV2GF(n e — 



\/2G F n e /2 + Acos20J + A 2 sin 2 26> 



(y/2G F n e /2 — A cos 29 



A 2 sin 2 26» 



(20) 



where a is +1 for the condition in Eq. 1)18(1 and —1 for 
the condition in Eq. (|19|) . Since in the Sun n e — n n > 
always, the above condition can be satisfied only for 
a = — 1 for 9 > 7r/4 and only for a = +1 for 9 < 7r/4. In 
other words, in the "dark side" 0,113 only the two light 
levels can intersect, while in the "light side", as preferred 
by the data, only the two heavy levels can. 

With this in mind, let us solve Eq. (|20|l . We find, in 
addition to the trivial vacuum solution, n e = n n — 0, a 
solution 



Am 2 
~2E^ 



= V2GF{n e -n n )i 



i 2 cos 2 29 — (n e 



(21) 



This is the desired general condition for having resonant 
antincutrino production in the solar magnetic field. Once 
again, for 9 < 7r/4, it solves Eq. I|18[l . while in the oppo- 
site case it solves Eq. (|19|) . (It is also worth reminding 
here that (|21|l has been derived assuming the neutrino is 
a Majorana particle.) 

We indeed see that Eq. (|21|l and the conventionally 
used Eq. @ agree only for 9 = and are substantially 
different for large 9. In particular, as 9 is increased be- 
yond 



| cos26> crit | ~ (1 - n n /n e ), 



(22) 



the level crossing as given by Eq. (|21|) disappears. It can 
be said that there is a "topological" difference between 



the cases of small and large mixing and Eq. l]22[l gives 
the angle at which the topology changes. 

The ratio n n /n e varies with the distance to the center 
of the Sun. Correspondingly, the numerical value of 6 CT it 
for which the disappearance takes place varies with A 
(which controls the location in the Sun where the level 
crossing could have happened). In the central region, 
n n /n e ~ 1/2, while for the outer regions n n /n e ~ 1/6. 
The corresponding variation in 9 cr n is 



tan 9 C 



(0.09- 0.33). 



(23) 



Let us finally discuss when each of the level-crossings, 
(|18|1 and (|19fl . can affect the neutrino evolution. The 
answer depends on the probability of finding the initial 
electron neutrino in each of the two neutrino mass eigen- 
states. This probability, in turn, depends on the fla- 
vor mixing angle at the production point. If the mat- 
ter potential y/2Gpn e dominates over the vacuum term 
Am 2 j(2E v \ the mixing angle at production approaches 
7r/2 and the neutrino is produced almost entirely in the 
heavy neutrino mass eigenstate. In the opposite limit, 
the neutrino is produced in the superposition of both 
mass eigenstates, with the probabilities cos 2 9 and sin 2 9. 

The rough boundary between these two regimes for 
the ~10 MeV 8 B neutrinos lies near Am 2 ~ 10~ 4 eV 2 : 
for smaller Am 2 including the best-fit region of Am 2 ~ 
8x 10 -5 eV 2 only the condition in Eq. I|18|l matters, while 
in the opposite case both Eq. I|18|l and Eq. I|19|) do. For 
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Regions of sensitivity to magnetic moment 




10" 1 1 10 

tan 2 9 



FIG. 4: The regions of the oscillation parameter space where one may expect the electron antineutrino flux above the KamLAND 
bound Q (three different shadings correspond to three different normalizations of the magnetic field). A very optimistic value 
of the transition moment, /i — 1 x 10~ 11 hb, was taken. Also shown are the regions allowed by the analysis of the KamLAND 
data (shaded regions in the middle) and the region allowed by the solar data (unfilled contours). 



Am 2 > 10~ 4 eV 2 we can thus expect sizable antineutrino 
production both in the "light" and "dark" sides, while for 
Am 2 -C 10~ 4 eV 2 the effect can happen only in the light 
side. 



D. Large mixing and u e production: numerical 
treatment 

To confirm this analysis, we have carried out a nu- 
merical scan of the oscillation parameter space, in the 
range 0.5 x 1CT 5 < Am 2 < 3 x 1CT 5 eV 2 . The calcu- 
lation was carried out for the realistic numerical solar 
density and chemical composition profiles |55j |. For the 
magnetic field profile, we have taken the lowest toroidal 
Ohmic eigenmode, as described in [25| and in detail in 
[4fj| . While higher eigenmodes could also be present, the 
choice of the lowest mode is expected to capture the es- 
sential physics: the field is axisymmetric, smooth, and 
vanishes at the RZ/CZ boundary. The only free parame- 



ter describing the field in this case is the overall mode nor- 
malization. We have performed scans for three different 
(highly optimistic) normalization values of the magnetic 
field, -B m ax = 0.5, 2, 5 MG, also taking a highly optimistic 
value of the transition moment, /i = lx 10 _11 /xs. 

The results of the scan are shown in Fig. 0] The three 
colored regions on the left side (small mixing) indicate 
the regions where the P e flux for the chosen normalization 
values of the magnetic field would exceed the KamLAND 
bound. The corresponding regions in the "dark" side, 
9 > it/ 4, are also shown. All the regions agree well with 
the analytical analysis of the previous Subsection. 

The Figure also shows the allowed regions from the 
KamLAND analysis of their reactor antineutrino spec- 
trum 0| , as well as the region that gives a fit to the solar 
neutrino data (unfilled contours). The best-fit ("LMA- 
I") region lies at the intersection of the solar and Kam- 
LAND regions and, as can be clearly seen, has no overlap 
with the regions where the v e flux would be measurable. 

We note that a similar conclusion was reached in 1241 . 
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using, however, an analytical criterion that differs from 
Eq. I|22l) . The relatively large v e flux cited in that work 
was obtained by assuming the field strength of 50 MG 
at the neutrino production point, which was only later 
shown to be not allowed lial. 



IV. CONVECTIVE ZONE 

Compared to the case of the RZ, the case of the CZ re- 
quires a different treatment. The difference has to do 
with the turbulent nature of the CZ magnetic fields. 
Below, for completeness, we briefly review the relevant 
physics. 

The main idea is that random variations |63l l(>'l| of 
the magnetic field on scales equal to, or smaller than, 
the neutrino oscillation length lessen the suppression of 
the oscillations by a large diagonal mass splitting. While 
the neutrino follows slow variations with scales longer 
than the oscillation length A osc ~ 1/A adiabatically, for 
shorter variations the adiabaticity condition is violated. 
The net result of the small-scale "noise" is that the neu- 
trino state vector performs a "random walk" in the oscil- 
lation space, gradually receding from its original position. 

To be specific, consider a simple model in which the 
magnetic field changes on scales of A corr < A osc . The 
probability of spin flip upon traversing the first correla- 
tion cell, A corr , is (/xSAcorr) (simple vacuum oscillation 
on scales smaller than the oscillation length). Since the 
fields in different A CO rr intervals are uncorrelated, as the 
neutrino traverses many domains the probabilities, not 
amplitudes, add up. For solar neutrinos, the experimen- 
tal bound on v e flux implies P < 1. So long as this is 
satisfied, the probability of spin flip after traveling the 
distance L can be written as P ~ (/Z-BA corr ) 2 L/A corr . 

For slow variations with A corr 3> A OS c, as already 
mentioned, the evolution is adiabatic and the conver- 
sion probability is exponentially suppressed. One can 
also imagine a case when the field domains are longer 
than Aosc, but have sharp boundaries. In such a case, 
the conversion probability is P ~ sin 2 20L/X COII ~ 
(/i£?/A) 2 L/A corr We see that in any case the conver- 
sion is most efficient for A corr ~ A osc ~ 1/A, namely, 
-Popt ~ (l^B) LX osc . 

Let us apply this to solar neutrinos. Once again, 
we are interested in the high-energy 8 B neutrinos, as 
these are the neutrinos the KamLAND bound applies 
to. These neutrinos are produced almost exclusively in 
the heavy mass eigenstate V2 and then propagate adi- 
abatically through the RZ. In the CZ, these neutrinos 
encounter turbulent magnetic fields and, if the transition 
moment is large, a small fraction of them gets converted 



into the light antineutrino mass eigenstate v%. The net 
result is the flux of neutrinos in the heavy eigenstate 
with possibly a small admixture of the antincutrinos in 
the light mass eigenstate. The light mass eigenstate, in 
turn, is detected as an electron antineutrino with proba- 
bility cos 2 8. The relevant field fluctuations are those on 
scales of A osc - 2 x 10 MeV/(8 x 10" 5 ) eV 2 - 3 x 10 2 
km. We will denote them by B\ osc . Overall, we get 

P{v e -> u e ) ~ cos 2 9(fxB Xo J 2 LX osc . (24) 

The oscillation parameters are known, cos 2 9 ~ 0.7, 
A - 10" 11 eV(10 eV/E v ). The extent of the region of 
the strong random field can be estimated to be some frac- 
tion of the CZ, say L ~ Q.IRq. Finally, the main issue 
is to estimate the size of the amplitude of the small-scale 
noise B\ asc . Unfortunately, the theory and simulation ef- 
forts are not yet at the level to give a robust prediction 
of -E>a osc (see for the state of the art on the simu- 
lation front). We will next consider two models for the 
magnetic field. 



A. Model of "uniform" turbulence 

A physically sensible estimate can be obtained by as- 
suming that (i) the field on the large scales is of the 
order of its turbulent equipartition value; and (ii) that 
the small-scale noise follows some kind of a scaling law 

B\ cx A Q , (25) 

ty pica l of turbulent systems. Such approach was adopted 
in |28ll29^ | and, in our opinion, is a clear improvement over 
naive models of delta-correlated noise. 

The largest scale of the turbulence, at which energy 
pumping takes place, is estimated to be of order the 
typical pressure scale height in the CZ. In the inner 
part of the CZ, the pressure scale height is of order 
imax ~ 5 x 10 4 km 44]. This number can be estimated 
directly from the solar model [55j, and is also consistent 
with the value of the mixing length 8 x 10 4 km at the 
base of the CZ given in [65j. Although the corresponding 
large-scale eddies reside in the interior of the CZ, they 
manifest themselves on the solar surface as supergranules 
(I ~ 15 x 10 4 km). These should not be confused with 
surface granules (I ~ 10 3 km) whose size is dictated by 
the scale height near the surface. 

The turbulent equipartition value for the large scale 
turbulent fields near the bottom of the CZ is estimated 
as 5 Lmax ~ i o 1 /6 i V3 r -2/ 3 _ io kG H. 

Putting the numbers together, we find 
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I" W B WW A osc \ 2a *( L \f E v \ f8xl0- 5 cV 2 \ fcos 2 
P(v e -ti/ e )~10 11-11 



10- n Ms/ V 10kG / V L max/ V 105 kn V V 10 MeVy V Am 2 /V 0.7 )' 

, (261 



This means that, for this field model, even with the most 
optimistic values of the neutrino magnetic moment Kam- 
LAND will not have a chance to see antineutrinos from 
the CZ spin flip unless its sensitivity is improved by as 
least one order of magnitude. 



B. Model of isolated flux tubes 

The estimate just given uses a simple model for the 
turbulent magnetic field. It is possible that the field in 
the CZ instead has a "fibril" nature ^4|, i- e -; * s expelled 
by the turbulence and combines in isolated flux tubes. It 
has been argued that the total energy of the CZ (thermal 
+ gravitational + magnetic) is reduced by the fibril state 
of the magnetic field by avoiding the magnetic inhibition 
of convection j6(|. 

Let us then consider a model in which the neutrinos 
travel through spatially separated flux bundles. The con- 
version probability in this case depends on the transverse 
size of a bundle. Given that a sunspot typically contains 
10 20 Mx [43L El and assuming superequipartition field 
strength inside the bundle of ~ 100 kG, one finds the 
transverse size of the bundle that is actually close to op- 
timal, ~ 300 km. The thickness of the flux rope comes 
out to be close to the oscillation length for the 8 B neu- 
trinos, a remarkable coincidence indeed! Since the total 
toroidal flux in the CZ is estimated to be of the order of 
10 24 Mx |43j, a given neutrino encounters only several of 
such flux bundles. The resulting conversion probability in 
the most optimistic case is then close to the KamLAND 
sensitivity bound, as can be easily checked. 



V. CONCLUSIONS 



to have large v e — > conversion even with magnetic 
fields an order of magnitude below the bound, B ~ 0.5 
MG [25|]. Moreover, as seen in Fig. 2] even a relatively mi- 
nor shift of the best-fit point, from Am 2 ~ 0.8 x 10~ 4 eV 2 
and ~ 32° to, e.g., Am 2 ~ 2 x 10~ 4 eV 2 and 9 ~ 24°, 
would have made the resonant antineutrino production 
possible. The field of solar neutrinos has truly entered 
the precision measurement stage! 

One may wonder if our analysis is completely self- 
consistent. After all, the oscillation parameters we use 
are found assuming only the flavor oscillations and ne- 
glecting possible magnetic v «-> v transitions. This, how- 
ever, is entirely justified, because the KamLAND bound 
on the P e flux ensures that the magnetic transitions can 
be present only at very small levels and thus cannot affect 
the fit to the oscillation parameters. 

In the convective zone, the solar neutrino propagate 
through turbulent magnetic fields. Our estimate of the 
probability of the spin-flavor flip in this case gives a rate 
that, at present, is likewise not observable. This conclu- 
sion disagrees with earlier analyses. The source of the 
disagreement is not in the treatment of the neutrino evo- 
lution, but in the treatment of the solar magnetic fields, 
which we believe were significantly overestimated in the 
earlier studies. 

Overall, our conclusion is that, despite achieving a 
clearly impressive level of sensitivity to P e from the Sun, 
KamLAND has not yet reached the level where it would 
be placing a meaningful constraint on the neutrino transi- 
tion moment and/or solar magnetic fields. On the other 
hand, KamLAND has reached the point where it may 
begin probing the "optimistic" scenario, \i v ~ 10~ n /j.B- 
Clearly, further sensitivity improvements through accu- 
mulation of additional statistics is very desirable. 



We have shown that the large measured value of the 
flavor mixing angle qualitatively changes the treatment of 
the neutrino interaction with the magnetic field in the ra- 
diative zone of the Sun. The physics of the effect is most 
transparent in the basis of the mass eigenstates, rather 
than in the flavor basis. We have shown that in the mass 
basis, for the allowed values of the oscillation parameters, 
there is no neutrino-antineutrino level crossing and, as a 
result, no measurable antineutrino production in the RZ. 
A correct level crossing condition, Eq. I|21|) . was derived, 
and a detailed numerical scan carried out to confirm this 
conclusion. 

We emphasize that the recent determination of the 
neutrino oscillation parameters is crucial for reaching this 
conclusion: for small flavor mixing angle, it was possible 
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APPENDIX A: SOME TECHNICAL DETAILS ON 
THE CONVECTIVE ZONE 

As described in Sect. IIVI simple physical arguments 
give the following estimates for the probability of spin- 
flavor flip in random magnetic fields: 

!(^A corr ) 2 L/A 
corr i Acorr ^ A 

corr; ^corr 
> A osc , sharp edge(Al) 
exp. suppressed, A cor r ^ A osc , smooth edge 

More generally, the probability of conversion is propor- 
tional to a Fourier component of the field correlation 
function 

/>oo 

P~H 2 L dl'(B(0)B(l'))cos(2Al'). (A2) 
Jo 

If (B(Q)B(l')) is a slowly decaying function (compared 
to 1/A) with a vanishing slope for small I', indicat- 
ing no small-scale noise, for example (B(0)B(l')) cx 
exp[— (r) 2 /A 2 orr ], the resulting P is exponentially small, 
as expected for adiabatic evolution. 

The quantity of interest in the isotropic turbulence 
model is B Xosc ~ £ Lmax (A osc /L max ) a . Refs. [Hill Per- 
formed an estimate based on a simple Kolmogorov turbu- 
lence. Within this model, the effects of the magnetic field 



on the character of the turbulence are neglected and the 
simple Kolmogorov scaling law for velocities is applied, 
v\ ~ (eA) 1 / 3 , where e is the energy dissipation rate per 
unit mass. The field was assumed to scale with A the 
same way as velocities, i.e., B\ cx A 1 / 3 . 

More generally, if the magnetic energy becomes com- 
parable to the kinetic energy of the fluid, the character 
of the turbulence changes (see, e.g., j65|, pp. 138-142). 
The rate of energy transfer to small scales is decided 
by the interaction of Alfven waves in plasma. Recent 
investigations argue that in the presence of mag- 
netic fields the turbulence becomes anisotropic in k space 
(see Appendix A, for a clear summary). While 

a detailed discussion of the latest models of magneto- 
hydrodynamics turbulence is clearly beyond the scope 
of this work, the preceding example show that the ex- 
ponent in the scaling law l|25|l is not known precisely. 
The correlation function can then be parameterized as 
(B(0)B(V)) cx (l'/L max ) a K a {l'/L max ), yielding upon in- 
tegration 

P ~ M 2 #2 max (A osc /L max ) 2Q L/ A, (A3) 
as expected on physical grounds. 
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While, by the late 1990's, the lack of time variations in 
the Super-Kamiokande data gave strong evidence against 
spin precession in the solar convective zone, spin preces- 
sion in the radiative zone continued to give a good fit to 
all solar data [25|. 

See also [kJ where some of the same arguments are men- 
tioned. 

Since the atmospheric neutrino measurements indicate 
large mixing in the v^ — v T sector, solar neutrinos actually 
oscillate between and a certain superposition state of 
and v T . To stay with the tradition, we denote the 
latter state by i/ M (and the corresponding superposition 
antineutrino state by 

The conversion probability in the resonance region 
depends on whether the cancellation persists "long 
enough", i.e., whether the oscillation length due to d fits 
inside the interval 81. This condition gives nothing but 
the usual adiabaticity parameter of the resonance |57fl. 
An extremely illuminating analysis of this in perturba- 
tion theory can be found in the classical work on the 
subject |60|. There, it was observed that the perturba- 
tive amplitudes for v e — > — > and z/ e — > — > 
have opposite signs and hence tend to cancel. The can- 
cellation is incomplete if one of the channels is enhanced 
with respect to the other one by matter effects. 
This was how the solution described in |25| was con- 
structed. Notice that there it was essential that the mix- 
ing angle was small. 



